Journal of Graph Algorithms and Applications 

http : / /www . cs . brown . edu/ publications/ jgaa/ 
vol. 3, no. 3, pp. 1-27 (1999) 



Subgraph Isomorphism in Planar Graphs and 
Related Problems 

David Eppstein 

Department of Information and Computer Science 
University of California, Irvine 
http : //www . ics .uci . edu/^eppstein/ 
eppsteinSics . uci . edu 

Abstract 

We solve the subgraph isomorphism problem in planar graphs in linear 
time, for any pattern of constant size. Our results are based on a tech- 
nique of partitioning the planar graph into pieces of small tree- width, and 
applying dynamic programming within each piece. The same methods 
can be used to solve other planar graph problems including connectivity, 
diameter, girth, induced subgraph isomorphism, and shortest paths. 
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1 Introduction 

Subgraph isomorphism is an important and very general form of exact pattern 
matching. Subgraph isomorphism is a common generalization of many impor- 
tant graph problems including finding Hamiltonian paths, cliques, matchings, 
girth, and shortest paths. Variations of subgraph isomorphism have also been 
used to model such varied practical problems as molecular structure compar- 
ison fl, integrated circuit testing |l0[] , microprogrammed controller optimiza- 
tion , prior-art avoidance in genetic evolution of circuits |^ , analysis of Chi- 
nese ideographs |2^ , robot motion planning , semantic network retrieval |^ , 
and polyhedral object recognition [Q. 

In the subgraph isomorphism problem, given a "text" G and a "pattern" 
H, one must either detect an occurrence of as a subgraph of G, or list all 
occurrences. For certain choices of G and H there can be exponentially many 
occurrences, so listing all occurrences can not be solved in subexponential time. 
Further, the decision problem is NP-complete. However for any fixed pattern 
H with i vertices, both the enumeration and decision problems can easily be 
solved in polynomial 0{n^) time, and for some patterns an even better bound 
might be possible. Thus one is led to the problem of determining the algorithmic 
complexity of subgraph isomorphism for a fixed pattern. 

Here we consider the special case in which G (and therefore H) are planar 
graphs, a restriction naturally occurring in many applications. We show that 
for any fixed pattern, planar subgraph isomorphism can be solved in linear 
time. Our results extend to some other problems including vertex connectivity, 
induced subgraph isomorphism and shortest paths. 

Our algorithm uses a graph decomposition method similar to one used by 
Baker to approximate various NP-complete problems on planar graphs. Her 
method involves removing vertices from the graph leaving a disjoint collection 
of subgraphs of small tree-width; in contrast we find a collection of non-disjoint 
subgraphs of small tree- width covering the neighborhood of every vertex. 

We assume throughout that all planar graphs are simple, so that the number 
of edges is at most 0{n); this simplifies our time bounds as we need not include 
the dependence on this number. The only problems for which this assumption 
makes a difference are induced subgraph isomorphism, ^.-clustering, and edge 
connectivity; for those, one can assume without loss of generality that the graph 
has bounded edge multiplicity, so again m = 0{n). 



2 New Results 

We prove the following results. The time dependence on H is omitted from 
these bounds. In general it is exponential (necessarily so, unless P=NP, since 
planar subgraph isomorphism is NP-complete) but see Theorem ^ for situations 
in which it can be improved. 

• We can test whether any fixed pattern H \sa, subgraph of a planar graph G, 
or count the number of occurrences of as a subgraph of G, in time 0{n). 
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• If connected pattern H has k occurrences as a subgraph of a planar graph 
G, we can hst all occurrences in time 0{n + k). If ff is 3-connected, then 
k = 0{n) and we can list all occurrences in time 0(n). 

• We can count the number of induced subgraphs of a planar graph G iso- 
morphic to any fixed connected pattern H in time 0(n), and if there are 
k occurrences we can list them in time 0{n + k). 

• For any planar graph G for which we know a constant bound on the 
diameter, we can compute the exact diameter in time 0{n). 

• For any constant h we can solve the /i-clustering and connected /i-cluster- 
ing problems in planar graphs in time 0{n). 

• For any planar graph G for which we know a constant bound on the girth, 
we can compute the exact girth in time 0{n). The same bound holds if 
instead of girth we ask for the shortest separating cycle or for the shortest 
nonfacial cycle in a given plane embedding of the graph. 

• For any planar graph G, we can compute the vertex connectivity and 
edge connectivity of G in time 0{n). (For planar multigraphs, we can test 
fc-edge-connectivity for any fixed k in time 0{n).) 

• For any planar graph G and any constant £, we construct in time 0{n) a 
linear-space routing data structure which can test for any pair of vertices 
whether their distance is at most £, and if so find a shortest path between 
them, in time 0{\ogn). 



3 Related Work 

For general subgraph isomorphism, nothing better than the naive exponential 
0(n'^') bound is known. Plehn and Voigt give an algorithm for subgraph 

isomorphism which in planar graphs takes time liJjC'd^Dn'^'v'l^^ (since im- 
proved by Alon et al. to 2<^(l^l)n°(Vl^)), but this is still much larger than 
the linear bound we achieve. 

Several papers have studied planar subgraph isomorphism with restricted 
patterns. It has long been known that if the pattern H is either K^^ or K^, then 
there can be at most 0{n) instances of as a subgraph of a planar graph G, 
and that these instances can be listed in linear time |2^, a fact which 
has been used in algorithms to test connectivity [p5[ , to approximate maximum 
independent sets [Q, and to test inscribability | [14| . Linear time and instance 
bounds for i^a and K^i can be shown to follow solely from the sparsity properties 
of planar graphs [|2[ and similar methods also generalize to problems of 
finding K2^2 and other complete bipartite subgraphs [^2|, |l^. Richards |Q gives 
O(nlogn) algorithms for finding G5 and Cq subgraphs in planar graphs, and 
leaves open the question for larger cycle lengths; Alon et al. Q gave 0{n log n) 
deterministic and 0{n) randomized algorithms for larger cycles. In ||l6|, we 
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showed how to list all cycles of a given fixed length in outerplanar graphs, in 
linear time (see also |3^, ^ for similar variants of outerplanar subgraph 
isomorphism) . We used our outerplanar cycle result to find any wheel of a given 
fixed size in planar graphs, in linear time. Itai and Rodeh |2|] discuss the 
problem of finding the girth of a general graph, or equivalently that of finding 
short cycles. The problem of finding cycles in planar graphs was discussed above. 
Fellows and Langston |Q discuss the related problem of finding a path or cycle 
longer than some given length in a general graph, which they solve in linear 
time for a given fixed length bound. The planar dual to the shortest separating 
cycle problem has been related by Bayer and Eisenbud |^ to the Clifford index 
of certain algebraic curves. Our results here generalize and unify this collection 
of previously isolated results, and also give improved dependence on the pattern 
size in certain cases. 

Recently we were able to characterize the graphs that can occur at most 0{n) 
times as a subgraph isomorph in an n- vertex planar graph: they are exactly the 
3-connected planar graphs ||l^ . However our proof does not lead to an efficient 
algorithm for 3-connected planar subgraph isomorphism. In this paper we use 
different techniques which do not depend on high-order connectivity. 

Laumond [ ^5[ gave a linear time algorithm for finding the vertex connec- 
tivity of maximal planar graphs. Eppstein et al. [ p^ give an 0(n) time algo- 
rithm for testing fc-edge-connectivity for fc < 4 and A;-vertex-connectivity for 
A; < 3. For general graphs, testing /c-edge-connectivity for fixed k takes time 
0(m -|- nlogn) [ p5[ . 4- vertex-connectivity in general graphs can be tested in 
time 0{na{n) -I- m) | p9| . However planar graphs can be as much as 5- vertex- 
connected, and nothing even close to linear was known for testing planar 5- 
connectivity. 

Our shortest path data structure combines our methods of bounded tree- 
width decomposition with a separator-based divide and conquer technique due 
to Frederickson ||2l| . Obviously all pairs shortest paths can be computed in time 
0{nm) after which the queries we describe can be answered in time 0(1], but 
some faster algorithms are known for approximate planar shortest paths ||2^, 

Our data structure answers shortest path queries exactly, in less prepro- 
cessing time than the other known results, but can only find paths of constant 
length. 

A final note of caution is in order. One should not be confused by the 
superficial similarity between the subgraph isomorphism problems posed here 
and the graph minor problems studied extensively by Robertson, Seymour, and 
others One can recognize path subgraphs by minor testing, but such tricks 
do not work for most other subgraph isomorphism problems. The absence of 
a fixed minor imposes severe structural constraints on a graph, whereas this is 
much less the case when a fixed subgraph is not present. Although minor testing 
can be done in time polynomial in the text graph size, the constant factors are 
typically much higher than those for our subgraph isomorphism algorithm. 
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Figure 1: Tree decomposition of a planar graph. 



4 Bounded Tree- Width Subgraph Isomorphism 

As a subroutine, we need to perform subgraph isomorphism testing in graphs 
of bounded tree-width. This can be done by a standard dynamic programming 
technique |9[ Q . The exact statement of the problem we solve is complicated by 
the requirement that we count or list each subgraph isomorph exactly once. For 
simplicity, we state the bounds for this problem with one parameter measuring 
both the tree-width of the text and the size of the pattern. 

Definition 1 A tree decomposition of a graph G consists of a tree T , in which 
each node N G T has a label L{N) C V{G), such that the set of tree nodes 
whose labels contain any particular vertex of G forms a contiguous subtree of T, 
and such that any edge of G connects two vertices belonging to the same label 
L{N) for at least one node N of T . The width of the tree decomposition is 
one less than the size of the largest label set in T . The tree-width of G is the 
minimum width of any tree decomposition of G. 

We can assume without loss of generality (by splitting high-degree nodes 
into multiple nodes with the same label) that each node in T has at most 
three neighbors and that there are 0{n) nodes in the tree. We will assign 
our tree decomposition an arbitrary orientation, by rooting it at one of its 
leaves, so that T becomes a binary tree. Figure ^ shows a planar graph, with 
a tree decomposition of width three. In fact, the graph shown has no tree 
decomposition of width two, so its tree-width is three. 

Define the subtree rooted at a node N to consist of N and all its descendants. 
Each such subtree is associated with an induced subgraph of G, having vertices 
contained in labels of nodes in the subtree. 
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Figure 2: Partial isomorph of a pentagon in the induced subgraph associated 
with node {F, H, K, A/}, and corresponding partial isomorph boundary mapping 
the pentagon to G"^. 

Lemma 1 The subtree rooted at N provides a tree decomposition of the asso- 
ciated induced subgraph of G. 

Proof: The only property of a tree decomposition that does not follow imme- 
diately is the requirement that each edge connect two vertices contained in the 
label of some node. Since this is true of G and T, any induced subgraph edge 
{u,v) must have {u,v} C L{N') for some N' , but N' may not be a descendant 
of N . However, if not, u belongs to both L{N') and (by assumption) L{N") 
where N" is a descendant of N . Therefore, by contiguity, u G L{N), and sim- 
ilarly V S L{N), so in this case {u,v) still both belong to the label of at least 
one node in the subtree. □ 

Lemma 2 Assume we are given graph G with n vertices along with a tree de- 
composition T of G with width w. Let S be a subset of the vertices of G, and 
let H be a fixed graph with at most w vertices. Then in time 2^^'^^°^^^n we 
can count all isomorphs of H in G that include some vertex in S. We can list 
all such isomorphs in time 2^^'^^°^'^^n-\-0{kw), where k denotes the number of 
isomorphs and the term kw represents the total output size. 

Proof: We perform dynamic programming in tree T. Let a partial isomorph 
at a node N of the tree be an isomorphism between an induced subgraph H' 
of the pattern H and the induced subgraph of G associated with the subtree 
rooted at N. 

We let G'j^ be formed by adding two additional vertices x, y to the subgraph 
of G induced by vertex set L{N). We connect each of the two additional vertices 



D. Eppstein, Planar Subgraph Isomorphism, JGAA, 3(3) 1-27 (1999) 7 



to all vertices in L{n), and each of the two additional vertices also is given a 
self-loop. Then from any partial isomorph at N we can derive a graph homo- 
morphism from all of H to G', which is one-to-one on vertices in L{N), maps the 
rest of H' to x, and maps H — H' to y. Let a partial isomorph boundary be such 
a map; Figure ^ illustrates a partial isomorph and the corresponding boundary. 
Since a partial isomorph boundary consists of a map from a set of at most w 
objects to a set of at most w + 3 objects, there are at most w'^'^^ = 2'^("''°s"') 
possible partial isomorph boundaries for a given node. 

Suppose that node N has children iVi and N2 We say that two partial 
isomorph boundaries B : H 1-^ G'j^ and Bi : H 1-^ G'^j are consistent if the 
following conditions all hold: 

• For each vertex w e iJ, if B{v) G L{Ni) or Bi{v) e L{N), then B{v) = 
B,iv). 

• For each vertex v e H, if B{v) ^ x then Bi{v) e L{N) U {y}. 

• At least one vertex v e H has Bi{v) ^ L{N) U {y}. 

We say that two partial isomorph boundaries Bi : H t-^ G'j^ and B2 : H G'^ 
form a compatible triple with B if the following conditions both hold: 

• Bi and B2 are both consistent with B. 

• For each v with B{v) — x, exactly one of Bi{v) and B2{v) is equal to y. 

For each partial isomorph boundary B : H i—^ G'^, let Xi{B) be the number 
of partial isomorphs which give rise to that boundary, and include a vertex of S. 
Let X2{B) be the number of partial isomorphs which give rise to that boundary, 
and do not include a vertex of S. These values can be computed in a bottom-up 
fashion as follows: 

• If there is no v for which B{v) = x, then all partial isomorphs having 
boundary B involve only vertices in L(N), and can be enumerated by 
brute force in time 20(«'i°g™). 

• Otherwise, we initialize Xi{B) and X2{B) to zero. Then, for each partial 
boundary Bi that is consistent with B, and such that there is no v with 
B{v) = X and Bi{v) = y, we increment Xi{B) by Xi{Bi) and increment 
X2{B) by X2{Bi). Finally, for each compatible triple B, Bi, B2 we incre- 
ment Xi{B) by Xi(Bi) • Xi{B2) + Xi{Bi) ■ ^2(52) + ^2(^1) • ^1(^2) 
and increment X2(B) by X2(-Bi) • X2(B2)- 

The total time for testing all triples for compatibility and performing the 
above computation is 0(w;3('"+3)-Hi = 2'^^'"^°s^)^ 

At the root node of the tree, we compute the number of isomorphs involving 
iS* simply by summing the values X2{B) over all partial isomorph boundaries 
for which B{v) 7^ y for all v. □ 



D. Eppstein, Planar Subgraph Isomorphism, JGAA, 3(3) 1-27 (1999) 8 



Lemma 3 Assume we are given graph G with n vertices along with a tree de- 
composition T of G with width w. Let S be a subset of the vertices of G, and 
let H be a fixed graph with at most w vertices. Then we can list all isomorphs 
of H in G that include some vertex in S in time 2'-'("' '°s™)ri + 0{kw), where 
k denotes the number of isomorphs and the term kw represents the total output 
size. 

Proof: We first follow the above dynamic programming procedure, to compute 
the values Xi and X2 for each partial isomorph boundary. We then compute 
top-down in the tree the set of pairs {B,X) where i? is a partial isomorph 
boundary and X is either Xi or X2, such that the value X{B) contributes to 
the final count of subgraph isomorphs. These pairs can be identified as the ones 
such that X{B) was included in the computation of some pair higher in the tree 
that has been previously identified as contributing to the total, and that caused 
a nonzero increment in this computation. Finally, we compute bottom-up again, 
listing for each contributing pair [B, X) the partial subgraph isomorphs counted 
in the value X{B). This step can be performed by mimicking the initial com- 
putation of X{B) described in the previous lemma, restricted to the boundaries 
known to contribute to the overall total, replacing each increment by a concate- 
nation of lists, and replacing each multiplication with the construction of partial 
isomorphs from a Cartesian product of two previously-computed lists. 

The number of steps for this computation is proportional to the number of 
steps in the previous algorithm, together with the added time for each combi- 
nation of a pair of partial isomorphs. Each such combination can be charged 
to a subgraph isomorph included in the output, and each output isomorph is 
formed by a binary tree of combinations that takes 0{w) time to perform, so 
the total added time is 0{kw). □ 

The same dynamic programming techniques also lead to similar results for 
counting or listing induced subgraphs isomorphic to H. To do this, we need only 
modify the algorithms above to restrict attention to partial isomorph boundaries 
B : H ^ G'j^ in which all edges between vertices of L{N) are covered by the 
image of some edge in H. 

5 Neighborhood Covers 

We have seen above that we can perform subgraph isomorphism quickly in 
graphs of bounded tree-width. The connection with planar graphs is the follow- 
ing: 

Lemma 4 (Baker [^) Let planar graph G have a rooted spanning tree T in 
which the longest path has length i. Then a tree decomposition of G with width 
at most 3£ can be found in time 0{in). 

Proof: Without loss of generality (by adding edges if necessary) we can assume 
G is embedded in the plane with all faces triangles (including the outer face). 
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Figure 3: Triangulated planar graph, with depth two tree T rooted at A (shown 
by heavy solid lines), and tree decomposition with nodes corresponding to faces 
of the graph and edges complementary to T. 



Form a tree with one node per triangle, and an edge connecting any two nodes 
whenever the corresponding triangles share an edge that is not in T (Figure ^ . 
Label each node with the set of vertices on the paths connecting each corner 
of the triangle to the root of the tree. Then each edge's endpoints are part of 
some label set (namely, the sets of the two triangles containing the edge), and 
the labels containing any vertex form a contiguous subtree (namely, the path 
of triangles connecting the two triangles containing the edge from the vertex 
to its parent, and any other triangles enclosed in the embedding by this path). 
Therefore, this gives us a tree decomposition of G. The number of nodes in any 
label set is at most 3^+1, so the width of the decomposition is at most 3£. □ 

In particular, any planar graph with diameter D has tree-width 0{D). 

If an isomorph of a connected pattern H uses vertex v in G, it is contained 
in the portion of G within distance \H\ of v. By Lemma ^ this \H\-neighborhood 
of V has tree- width at most 3\H\. Therefore we can cover G by the collection 
of all such neighborhoods, and use Lemma ^ to find the copies of H within 
each neighborhood. However such a cover is not efficient: the total size of all 
subgraphs is O(n^), so this would give us a subgraph isomorphism algorithm 
with quadratic runtime. We speed this up to linear by using more efficient 
covers. 

Awerbuch et al. |^ have introduced the very similar concept of a neigh- 
borhood cover, which is a covering of a graph by a collection of subgraphs, with 
the properties that the neighborhood of every vertex is contained in some sub- 
graph, and that every subgraph has small diameter. They showed that for any 
(possibly nonplanar) graph, and any given value w, there is a w-neighborhood 
cover in which the diameter of each subgraph is 0{w\ogn), and in which the 
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total size of all subgraphs is 0{m\ogn); such a cover can be computed in time 
0(m log n + nlog^ n) [Q. Because of Lemma ^, such a neighborhood cover is 
also almost exactly what we want to speed up our subgraph isomorphism algo- 
rithm. However there are two problems. First, the size and construction time 
of neighborhood covers are higher than we want (albeit only by polylogarithmic 
factors). Second, and more importantly, the diameter of each subgraph is loga- 
rithmic, so we are unable to use dynamic programming directly in the subgraphs 
of the cover. We would instead be forced to use some additional techniques such 
as separator-based divide and conquer, introducing more unwanted logarithmic 
factors. 

Instead, we use a technique similar to that of Baker ||^ to form a cover that 
has the properties we want directly: any connected lu-vertex subgraph of G is 
included in some member of the cover, and each vertex of G is included in few 
members of the cover (so the total size of the cover is 0{n)). Unlike the tech- 
niques cited above, the diameter of the subgraphs will not be bounded, however 
we will still be able to use Lemma ^ on an auxiliary graph to show that each 
covering subgraph has tree- width 0(w). Because of the exponential dependence 
of our overall algorithms on the tree-width of the covering subgraphs, we con- 
centrate our efforts on reducing this width as much as possible, at the expense 
of increasing the total size of the cover by an 0{w) factor over the minimum 
possible. 

Lemma 5 Let G be a planar graph, and w be a given integer parameter. Then 
we can find a collection of subgraphs Gi and a partition of the vertices of G into 
subsets Si with the following properties: 

• Every vertex of G is included in at most w subgraphs Gi . 

• We can find a tree decomposition of each subgraph Gi with tree-width at 
most 3w — 1. 

• If H is a connected w-vertex subgraph of G, and i is the smallest value 
for which H C] Si is nonempty, then H is a subgraph of Gi but is not a 
subgraph of any Gj with j > i. 

• The total time for performing the partition and computing the tree decom- 
positions is 0{w^n). 

Proof: We choose an arbitrary starting vertex uq, and let Si consist of the 
vertices at distance i from vq. We then let Gi be the graph induced by the 
vertex set ljj=r~^ shown in Figure^. Clearly, the sets Si form a partition 

of the vertices of G, and each vertex is in at most w subgraphs Gi. 

Then for i = 0, Gi consists of the vertices at distance at most w — 1 from 
vq, so by applying Lemma ^ to its breadth first spanning tree we can find a 
tree decomposition with width at most 3(w — 1). To show that each Gi with 
i > has low tree-width, form an auxiliary graph G[ from G by collapsing into 
a single supervertex all the vertices at distance less than i from wq, and deleting 
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Figure 4: Planar graph with breadth first spanning tree (heavy edges), partition 
into layers Si, subgraph Gi (for w = 3, i — 5), and minor G[ (with large 
supervertex and contracted breadth first spanning tree). 



all the vertices with distance at least i + w. is a minor of the planar graph G 
and is therefore also planar. Then similarly collapsing a breadth first spanning 
tree of G gives a spanning tree of with depth at most w, so G^ has a tree 
decomposition with width at most Sw, in which each node of the decomposition 
includes the collapsed supervertex in its label. Gi is formed by deleting this 
supervertex from G^, so we can form a tree decomposition of Gi with width at 
most 3w — 1 by removing the supervertex from the decomposition of G'j . 

Next, we need to show that any connected subgraph H of G with \H\ < w 
is contained in Gi, where i is the smallest value such that H O Si 7^ 0. But Gi 
is formed from G simply by removing the sets Sj where j<i or j>i + w. No 
Sj with j < i can contain a vertex of H, or else i would have been smaller. And 
no Sj with j > i + w can contain a vertex v of H, or else we could find a path 
of length at most i + w — 1 from vq to v by concatenating a path in H from 
some vertex Vi G Si H to v (which has length at most \H\ — 1 < w — 1) with 
the breadth first tree path from vq to Vi (which has length i), contradicting the 
placement of v in 5^-. Therefore, none of the vertices that were deleted from G 
can belong to H, so H remains a subgraph of Gi. 

Finally, the condition that H can not be a subgraph for Gj where j > i is 
clearly true, since no such Gj can include any vertex of Si. 

The time bound is dominated by the time to perform the tree decompositions 
on the graphs G-, which by Lemma § is ^ 0{wGi) = wO{J2 Gi) = w{wn). □ 
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6 The Subgraph Isomorphism Algorithm 

We first describe the result for tlie special case of connected patterns. 

Theorem 1 We can count the isomorphs or induced isomorphs of a given con- 
nected pattern H , having w vertices, in a planar text graph G with n vertices, 
in time 

20(wiogw)^^ If there are k such isomorphs we can list them all in time 
20i^^oE^)n^0iwk). 

Proof: The algorithm consists of the following steps: 

1. Apply the method of Lemma ^ to find a partition of the vertices into sets 
Si associated with graphs Gi having low width tree decompositions. 

2. For each i > 0, count or list the subgraph isomorphs of H in Gi that involve 
at least one vertex of Si, using the algorithm of Lemma ^ or Lemma |^ 
respectively. 

3. Sum all the counts or concatenate the lists, to get a count or list of the 
isomorphs in G. 

By Lemma m each isomorph of in G occurs in exactly one way as an isomorph 
in Hi that involves at least one vertex of Si, so the algorithm produces the 
correct total count or list. The time for the first step is 0{w^n), and the time 
for the last step is 0{n), both dominated by the time for the second step which 
is J]20("'i°s™)|Gj| +0(fc,w), where the 20("'i°8'"' factor arises by plugging the 
2iw — 1 treewidth bound of Lemma |^ into the analysis in Lemmas I and I 
This can be simplified to 20('"i°g«')(^ \G.,\) + 0(Y,hw) = 20('"'°s») • ©(wn) + 
O(fcw) = 2O(™>°s™)n + 0(M. □ 

The method so far requires that the pattern be connected. We now describe 
a general method for handling disconnected patterns. The technique will let 
us count the number of matching patterns, after which some sort of separator- 
based divide and conquer can likely be used to find an instance of a matching 
pattern, but we have been unable to extend this technique to the problem of 
listing all subgraph isomorphs of a disconnected pattern. 

Theorem 2 We can count the isomorphs of any (possibly disconnected) pattern 
H having at most w of vertices, in a planar text graph G with n vertices, in time 

20(tulogtu)^^ 

Proof: Let i^oiH) denote the number of isomorphs of H in G. Rather than 
counting the isomorphs of a single pattern, we count the isomorphs of all planar 
graphs having at most w vertices. There are only 2*^'™) such graphs [^, so 
this factor does not change the overall form of our time bound. We order these 
graphs by the number of connected components, so that when we are processing 
a particular graph H we can assume we already know the values of ^g{H') for 
every H' with fewer components. 

Our algorithm them performs the following steps on each graph H: 
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1. If is connected, compute ^q{H) using the algorithm of Theorem [l|. 

2. Otherwise, let H be the disjoint union of two subgraphs H' and H" , and 
let #g{H) = ^g{H') ■ ^g{H") - Y: h#GiH.,), where the sum is over all 
graphs Hi with fewer components than H , and ki denotes the number of 
different ways Hi can be formed as the union of Hi and H2. 

The product ^g{H') ■ ^g{H") counts the number of ways of mapping H 
into G such that both H' and H" are isomorphically mapped but their instances 
may overlap. The term ^^ki^GiHi) corrects for these overlaps by subtracting 
the number of overlapped occurrences of each possible type. 

The coefficients ki may be computed by brute force enumeration of all pos- 
sible ways of marking a vertex of Hi as coming from H' , H" , or both, combined 
with a planar graph isomorphism algorithm, in time 2'-'^''"K Therefore, the over- 
all time taken in the second step of the algorithm is 2'^^^\ independent of n, 
and the total time is dominated by the first step, in which we apply Theorem^ 
to 20("') connected graphs, taking time 2C('"i°g"')n. □ 



7 Further Improvements 

For certain patterns, such as the wheels, our results can be further improved 
to reduce the time dependence on \H\. Let diam(i7) denote the diameter of H 
(i.e., the longest distance between any two nodes), and let Kx{H) denote the 
maximum number of connected components that can be formed by removing 
at most X nodes from H . Note that if the diameter diam(iJ) is small, we can 
use that value instead of \H\ in our neighborhood cover of G, reducing the 
tree-width of the subgraphs Gi to 0(diam(i?)). 

Lemma 6 Let H he a given pattern graph, and N be a node of a tree decompo- 
sition of graph G. Then there are at most 2^ii(")i l+i) different 
possible partial isomorph boundaries B : H G'j^ (as defined in the proof of 
Lemma |^. 

Proof: The map B can be defined by specifying which (if any) vertex of H 
maps to each vertex in L{N) (using log(|i/| + 1) bits per vertex of L{N) to 
specify this information) and also specifying which of the remaining vertices in 
H map to the vertex x in and which ones map to the vertex y. However, it 
is not possible for the boundary to come from an actual subgraph isomorphism 
unless each connected component oi H \ B~^ {L{N)) is mapped consistently 
either to x or to y, since any path from a: to y must pass through a vertex of 
L{N). So, to finish specifying the boundary, we need only add this single bit of 
information per component oi H \ B~^{L{N)), and by definition there are at 
most iir|i(jv)|(^^) such components. □ 

As a consequence, the analysis in Lemma can be tightened to show that 
the dynamic program takes time 2'^^^oid.an.(H)^H)+dmm{H)iog\H\)^_ 
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Lemma 7 Suppose that planar graph H is Hamiltonian, or is 3- connected, or 
is connected and has hounded degree. Then for any set S of vertices of H , H — S 
has 0{\S\) connected components. 

Proof: For Hamiltonian graphs and bounded degree graphs this is straightfor- 
ward. For 3-connected graphs, assume without foss of generahty that no edge 
can be added to H connecting two vertices in S; then each component of H — S 
must occupy a distinct face in the planar embedding of S induced by the unique 
embedding of i?. □ 



Theorem 3 If a given pattern H is Hamiltonian, 3-connected, or connected of 
bounded degree, we can count the isomorphs of H in a planar text graph G with 
n vertices in time 20(diam(ff) log 

Proof: The proof consists simply of plugging the improved analysis of Lemma ^ 
into the algorithm of Theorem |l|. □ 

For instance we can count the isomorphs of a wheel Wk in a planar text 
graph G with n vertices, in time 0(nfc'^^'^^). In fact in this case it is not difficult 
to come up with an 0{nk^) algorithm: 

Theorem 4 We can count the isomorphs of any wheel Wk in a planar text 
graph G with n vertices in time Oink^). 

Proof: For each vertex v, we count the number of cycles of length k in the 
neighbors of v. The sum of the sizes of all neighborhoods in G is 0(n). Each 
neighborhood has treewidth at most 2 by Lemma ^ Any partial isomorphism 
of a fc-cycle in a node N of this decomposition can only consist of a single path 
of at most fc — 1 vertices, which starts and ends at some two of the at most 
three vertices in L{N) and may or may not involve the third vertex; therefore 
we need only keep track of 0{k) different partial isomorph boundaries at each 
node. A careful analysis of the steps in the algorithm of Lemma || then shows 
that the most expensive step (finding compatible triples) can be performed in 
time O(fc^) per node, giving an overall running time of 0(k^n). □ 



8 Variations and Applications 

We now describe briefly how to use our subgraph isomorphism algorithm to 
solve certain other related graph problems. For instance, we can find the girth 
(shortest cycle), or smallest nonfacial cycle (for some particular embedding) 
simply by searching for isomorphs of small cycles: 

Theorem 5 We can find the girth g of any planar graph G, or find a small- 
est nonfacial cycle G for an embedded planar graph, in time 2'^'^^'°s9)n or 
20(|c|iog|c|)^ respectively. 



D. Eppstein, Planar Subgraph Isomorphism, JGAA, 3(3) 1-27 (1999) 15 



Proof: We test for each integer i = 3, 4, 5 . . . whether there is a cycle (or 
nonfacial cycle) of length i. To test if there is a nonfacial cycle, we count the 
total number of cycles of length i in the graph and subtract the number of faces 
of length i. The total time for this procedure is J2z<g 2'="^'^°s'')n = 20(si°g9)n. 
For the nonfacial cycle problem, once the length of the cycle is known, we can 
find a single such cycle by performing our subgraph isomorph listing algorithm, 
stopping once k + 1 cycles are generated, where k = 0{n) is the number of faces 
of the given length. By radix sorting the list of cycles (in lexicographic order by 
their sequences of vertex indices) we can then test in linear time which of the 
generated cycles are nonfacial. □ 



Theorem 6 We can find a shortest separating cycle C in a planar graph, in 
time 20(tc|iog|c|)„^ 

Proof: We describe how to test for the existence of a separating cycle of length 
i; the shortest such cycle can then be found by a sequential search similar to 
the computation of the girth. 

We first modify the construction of the graphs Gi in Lemma |^, by including 
in Gi not only all the vertices in layers Si through Si+w-i but also a single 
supervertex for each connected component of the graph induced by the vertices 
in layers Sj, j > i + w, and a supervertex for the (single) connected component 
of the vertices in layers Sj, j < i. Then a cycle that uses only vertices in layers 
Si through S'i+uj-i (and does not use any of the supervertices) is separating in 
the modified Gi if and only if the corresponding cycle is separating in G. Note 
that the added supervertices only add one level to the breadth first search tree 
of Gi and hence the tree- width is still 0{w). 

Then we need merely modify the dynamic program of Lemma ||, to use a 
definition of a partial isomorph boundary that, in addition to the map B : 
H l—^ G'j^, specifies which of the remaining unmapped vertices of are in 
each of the two subsets of G separated by the cycle, enforcing the requirement 
that no vertex in one subset be adjacent to a vertex in the other subset. This 
modification multiplies the number of boundaries by 2^^'^\ but this increase is 
swamped by the 2'-'(™'°s™) terms from our previous analysis. □ 

We next consider the application of our techniques to certain graph clustering 
problems. 



Definition 2 (Keil and Brecht |30|) The /i-clustering problem is that of ap- 
proximating the maximum clique by finding a set of h vertices inducing as many 
edges as possible. The connected /i-clustering problem adds the restriction that 
the induced subgraph be connected. 

Keil and Brecht study these problems, and show that even though cliques 
are easy to find in planar graphs [40|, the connected ^.-clustering problem is NP- 
complete for planar graphs. See [32| for approximate /i-clustering algorithms in 
general graphs. 
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Figure 5: An embedded planar graph G, the vertex connectivity substitute 
graph G' (with edges drawn as heavy curves), and the edge connectivity substi- 
tute graph G". 

Theorem 7 For any h we can solve the planar h-clustering and connected h- 
clustering problems in time 2'-'^^^°^'^^n. 

Proof: We simply to test subgraph isomorphism for all possible planar graphs 
on h vertices, and return the subgraph isomorph with the most edges. □ 

We now describe two applications to connectivity that, unlike the previ- 
ous applications, are linear without an exponential dependence on a separate 
parameter. 

The vertex connectivity of G is the minimum number of vertices such that 
their deletion leaves a disconnected subgraph. Since every planar graph has a 
vertex of degree at most five, the vertex connectivity is at most five. We now 
use a method of Nishizcki (personal communication) to transform the vertex 
connectivity problem into one of finding short cycles, similar to those discussed 
at the start of this section. We choose some plane embedding of G, and construct 
a new graph G' having n + f vertices: the n original vertices of G and / new 
face-vertices corresponding to the faces of G. We place an edge in G' between an 
original vertex and a face- vertex whenever the corresponding vertex and face are 
incident in G. Then G' is a bipartite plane-embedded graph. This construction 
is illustrated in the center of Figure |[ 

Lemma 8 Any minimal set S d G of vertices the deletion of which would 
disconnect G corresponds to a cycle G <Z G' of the same original vertices and 
an equal number of face-vertices in G' , such that G'\C has at least two connected 
components each containing at least one original vertex. Conversely if G is any 
cycle such that G' \ G has at least two connected components each containing at 
least one original vertex, then the original vertices in G form a cutset in G. 
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Proof: Let A be a connected component oi G\ S, let Ga be formed from G 
by contracting A into a supervertex, and let 5" be the set of faces and vertices 
adjacent to the contracted supervertex. Then (since it is just the neighborhood 
of a vertex) 5" has the structure of a cycle in G", and separates A from G\{AUS}. 
If S is minimal, then it must consist of exactly the original vertices in S'. The 
converse is immediate, since no edge in the embedding of G can cross a face or 
vertex in C. □ 



Theorem 8 We can compute the vertex connectivity of a planar graph in 0{n) 
time. 

Proof: We can assume without loss of generality that G is two-connected, so 
the graph G' described in Lemma || has no multiple adjacencies. We form G" 
as above and find the shortest cycle in G" that separates two original vertices. 
As with the shortest separating cycle problem (Theorem ^) this can be done by 
a slight modification to our dynamic programming method that decorates the 
dynamic programming states with 0(1) bits of additional information regarding 
the separated vertices. Since any planar graph has a vertex of degree at most 
five by Euler's formula, the shortest cycle in G' must have length at most ten, 
so the algorithm takes time 0{n). □ 

The edge connectivity of a graph is similarly defined as the minimum number 
of edges the removal of which disconnects the graph. For simple graphs, this 
can again be at most five but for multigraphs it can be higher. 

Theorem 9 We can compute the edge connectivity of a simple planar graph in 
Gin) time. 

Proof: Assume without loss of generality that G is two-edge-connected. Em- 
bed G, and form a graph G" by subdividing each edge of G and connecting 
the resulting subdivision points to new vertices in each adjacent face; this con- 
struction is illustrated in the right of Figure Then G" is a planar graph 
with n original vertices, e edge-vertices on each edge of G, and / face-vertices 
in each face of G, so its total complexity is 0{n). By the assumption of two- 
edge-connectivity, G" is simple. One can use an argument similar to the one 
in Lemma |[ in which we delete a cutset from G, contract a connected com- 
ponent, and examine the neighborhood in G of the contracted supervertex, to 
show that G is fc-edge-connected iff there is no cycle of fewer than 2k edge- and 
face- vertices in G", which separates two original vertices. As before, the degree 
bound on a planar graph imposes a limit of ten on the length of the shortest 
such cycle, and as before this cycle can be found by a minor modification to our 
dynamic programming algorithm. □ 

The same methods extend easily enough to multigraphs, but now the edge 
connectivity can not be bounded a priori, so we need to include the connectivity 
in our time bound. 



D. Eppstein, Planar Subgraph Isomorphism, JGAA, 3(3) 1-27 (1999) 18 



Theorem 10 For any fixed k, we can test k- edge- connectivity of a planar multi- 
graph in time 2'^^^^°^^^n. 

Proof: Without loss of generality, the multiplicity of any edge is at most fc, as 
higher multiplicities can not improve the overall connectivity. After the edge 
subdivision step in the construction of G", the resulting graph is a simple planar 
graph with 0{kn) vertices, after which we can proceed as in the remainder of 
Theorem |. □ 



9 Shortest Path Data Structure 

We next describe a technique for finding shortest paths in planar graphs. Let a 
parameter ^ be given (typically, a fixed constant). We wish to test, for any two 
vertices u and w, whether there is a path from u to w of distance at most I, and 
if so return the shortest such path. 

Since we wish to use an amount of space independent of we need a variant 
of Lemma ^ in which the total size of the subgraphs is not so large. 

Lemma 9 Let G he a planar graph, and w be a given integer parameter. Then 
we can find a collection of subgraphs Gi and a partition of the vertices of G into 
subsets Si with the following properties: 

• Every vertex of G is included in at most two subgraphs Gi . 

• We can find a tree decomposition of each subgraph Gi with tree-width 
0{w). 

• Gi contains the £-neighborhood of every vertex in Si. 

• The total time for performing the partition and computing the tree decom- 
positions is 0{wn). 

Proof: As in Lemma ^, we compute the distances of each vertex from some 
arbitrary starting vertex vq. We then let Si consist of those vertices with dis- 
tance at least 2iw and at most (2i + 2)w — 1 from vq, and we let Gi be the 
graph induced by the set of vertices with distances at least (2i — l)w and at 
most (2i + 3)w — 1 from vq. The proof that these graphs have treewidth 0{w) 
and that each Gi contains the ^-neighborhood of Si is essentially the same as 
that of Lemma |^. □ 

As before, by introducing dummy nodes, we can assume without loss of 
generality that each node in the tree decomposition of Gi has at most three 
neighbors. We warm up with a data structure for our shortest path queries that 
uses more space than necessary, but for which queries are very fast. 

Theorem 11 Given a planar graph G, and any value £, we can in time and 
space 0{inlogn) build a data structure that can, given a query pair of vertices, 
either return the distance between the pair or determine that the distance is 
greater than i, in time 0{£) per query. 
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Proof: By performing the decomposition of Lemma Q we can assume without 
loss of generality that we have a tree decomposition T for G of width 0{£). 
As with any tree, we can find a node N the removal of which disconnects T 
into subtrees of size at most \T\/2. Our primary data structure consists of the 
distances d{x, z) from each vertex x G G to each vertex z G L{N), together with 
a recursively constructed data structure in each subtree. 

To answer a query pair x, y where the two vertices belong to different subtrees 
of T, we can simply try each of the 0{D) values d{x, z) +d{z, y) where z ranges 
over all the members of L{N). To answer a query where the two vertices belong 
to the same subtree, we can use the recursively defined structure in that subtree. 

It remains to show how we quickly determine which node N is eventually 
used to answer each query. To do this, define the level of a node to be the 
stage in the recursive subdivision process at which the node was chosen, and 
define the superior of a node N to be the node chosen at the next earlier level 
in the subtree containing N. The links from a node to its superior define a 
tree structure T' different from the original decomposition tree T. Further, 
define the home node of a vertex v to be the node with the earliest level with 
V e L{N). Note that, because of the requirements that the labels containing v 
are contiguous, the home node is uniquely defined. Then, the node to be used 
in answering a query pair x,y is simply the least common ancestor in T' of the 
home nodes of x and y. □ 

To return the actual shortest path, rather than simply the distance between 
a pair of nodes, we can store a single-source shortest path tree for each member 
of L{N), and return the path in the tree for the member of L(N) giving the 
smallest distance. 

We next show how to reduce the space to linear, at the expense of increasing 
the query time. 

Theorem 12 Given a planar graph G, and any value we can build a data 
structure of size 0{in) that can, given a query pair of vertices, either return the 
distance between the pair or determine that the distance is greater than £, in 
time 0{£^ \ogn) per query. 

Proof: As above, we can assume G has a tree decomposition of width 0(£), 
which we assume has the form of a rooted binary tree. Define levels in this tree 
and home nodes of vertices as above, except that we terminate the recursive 
subdivision process when we reach subtrees with fewer than i nodes (which we 
call small subtrees). If the node labels containing a vertex v belong only to 
nodes in a small subtree S, then v does not have a home node, instead we say 
that S is v's home subtree. 

Define a pair of nodes in T to be related if there is no node between them 
with an earlier level than both. Then, each node N is related to 0(1) nodes at 
each level: one node at an earlier level than N, and at most one node in each 
later level in each of the at most three subtrees formed by removing N from T. 
Therefore, there are 0{n/£) pairs of related nodes. 
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Our data structure then consists of the matrix of distances from vertices in 
L{Ni) to vertices in L{N2), for each pair N1,N2 of related nodes. The space 
for this data structure is 0{n£). It can either be built as a subset of the data 
structure of Theorem |l^, in time 0{ni\ogn), or bottom-up (using hierarchical 
clustering techniques of Frederickson to construct the level structure in 
T, and then computing each distance matrix from two previously-computed 
distance matrices in time 0{i^)) in total time 0{n£'^); we omit the details. 

To answer a query, we form chains of related pairs connecting the home 
nodes (or small subtrees) of the query vertices to their common ancestor in T'. 
The levels of the nodes in these two chains becomes earlier at each step towards 
the common ancestor, so the total number of pairs in the chain is O(logn). We 
then build a graph, in which we include an edge between each pair of vertices 
in the labels of a related pair of nodes, labeled with the distance stored in the 
matrix for that pair. We also include in that graph the edges of G belonging to 
the small subtrees containing the query vertices, if they belong to small subtrees. 
The query can then be answered by finding a shortest path in this graph, which 
has 0{£\ogn + P) vertices and 0{t^ logn) edges. □ 

The following theorem on computing diameter improves the naive O(n^) 
bound for all pairs shortest paths when the diameter is small. Note that diam- 
eter is not a subgraph isomorphism problem but it succumbs to similar tech- 
niques. 

Theorem 13 We can compute the diameter D of a planar graph G, in time 

Proof: We begin by performing a breadth first search from an arbitrary vertex. 
This will produce a tree of height at most D, so by Lemma ^ we can find a tree 
decomposition of width 0{D), which as usual we can assume has the form of a 
rooted binary tree. We first perform a bottom- up sweep of this tree to compute 
for every node N the distances between every pair of vertices in L{N), in the 
graph associated with the subtree rooted at N. These 0{D^) distances can 
be found by combining the distance matrices of the two children of N , in time 
0{D^), so this phase takes time 0{D^n). We then sweep the decomposition top 
down, computing for every node N the distances between every pair of vertices 
in L{N), in the whole graph G. The first pass correctly computed these distances 
at the root of the tree, and at any other node N the distances can be computed 
by combining the distance matrices of the parent of N (previously computed in 
the top-down sweep) and its two children (computed in the bottom- up sweep), 
again in time 0{D^) per node. 

We finally sweep through the tree decomposition bottom up, keeping at each 
node N a subset S of the vertices seen so far in the subtree rooted at N, together 
with the distances from each member of S to each member of L{N). When we 
process a node iV, we perform the following steps: 

1. Let the set S for node N consist of the union of the corresponding sets 5*1 
and 52 for its children A^i and N2, together with L{N). 
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2. Compute the distances from each member of S to each member of L{N), 
by combining the previously computed distances to L{Ni) or L{N2) with 
the distances within L{Ni) U L{N2). 

3. Compute the distance between each pair of nodes x, y where x ^ Si \ L(N) 
and y £ S2 \ L{N), by testing the distances through aU 0{D) possible 
intermediate nodes in L{N). 

4. Radix sort the members of S according to the lexicographic ordering of 
their 0(_D)-tuples of distances to L{N), and eliminate all but one member 
for each distinct tuple. 

The value returned as the diameter is then the maximum of the distances from 
S to L{N) computed in the second step, and the distances of pairs x, y comuted 
in the third step 

If any eliminated member x of a tuple belongs to a diametral pair x, y, where 
y is not in the subgraph associated with TV, then the uneliminated member x' 
with the same tuple would have the same distance to y, and would form another 
diametral pair. Therefore, the algorithm above will correctly find and report a 
diametral pair. 

The number of distinct 0(£')-tuples of integers in the range from to D 
is 2'^(^'°s^)j hence this gives a bound on the size of each set S. The time to 
compute distances between pairs x,y is 0{D) times the square of this quantity, 
which is still 20(^i°s^). □ 



10 Other Graph FamiUes 

Our results for planar graphs use the assumption of planarity in two ways. First, 
in the bound relating tree- width to diameter (Lemma ^), the proof is based on 
the existence of a planar embedding of the graph, and in fact there is no similar 
bound in general for nonplanar graphs; for instance the complete graph Kn has 
diameter one but tree- width n — 1. Second, in the cover of G by low-treewidth 
subgraphs described in Lemma ^, we needed the fact that planar graphs are 
closed under minors to show that the graph G'^ is planar, allowing us to apply 
Lemma ^ to it. 

This naturally raises the question, for which other minor-closed graph fami- 
lies can we prove a bound relating diameter to tree- width, similar to Lemma ^ 
Such a result would then let us apply our subgraph isomorphism techniques 
unchanged to any such families. In the conference version of this paper |l8| , we 
announced an exact characterization of these families, which are detailed in a 
separate journal paper ||l^ and which we now summarize: 

Definition 3 Define a family T of graphs to have the diameter-treewidth prop- 
erty if there is some function f{D) such that every graph in T with diameter at 
most D has treewidth at most f{D). 
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Figure 6: The graph on the left is an apex graph; the topmost vertex is one of 
the possible choices for its apex. The graph on the right is not an apex graph. 



Definition 4 An apex graph is a graph G such that for some vertex v (the 
apexj, G \ {v} is planar (Figure 



Theorem 14 (|15|) Let T be a minor- closed family of graphs. Then T has 
the diameter-treewidth property iff T does not contain all apex graphs. 



Corollary 1 Let T be a minor-closed family of graphs, such that some apex 
graph is not in T . Let H be a fixed graph in T . Then there is a linear time 
algorithm for testing whether H is a subgraph of a given graph G G JF. 

Note that the bound on tree-width from Theorem ^ is much higher than 
the linear bound in Lemma ^, so the dependence on of the time bound of 
the algorithm implied by Corollary ^ is much greater than in our planar graph 
algorithms. However, for certain important minor-closed graph families (such as 
bounded genus graphs) we were able to prove a better dependence of tree- width 
on diameter |l5[, leading to less impractical algorithms. 



11 Conclusions and Open Problems 

We have shown how to solve planar subgraph isomorphism for any pattern in 
time 0{n). We have also solved certain related problems in similar time bounds. 
A number of generalizations of the problem remain open: 

• We have shown that we can solve planar subgraph isomorphism even for 
disconnected patterns in time 0(n). Can we list all occurrences of a 
disconnected pattern in time 0{n-\- k)l 
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• Bui and Peck |Tl] describe an algorithm for finding the smallest set of 
edges partitioning a planar graph into two sets of vertices with specified 
sizes; if the edge set has bounded size their algorithm has cubic running 
time. Can we use our methods to find such a partition more quickly? 

• It seems possible that the recently discovered randomized coloring tech- 
nique of Alon et al. M can improve the dependence on the size of the 
pattern from 2'^'^"'^°'^^"' to 2^'^^\ but only for the decision problem of sub- 
graph isomorphism. Can we achieve similar improvements for the counting 
and listing versions of the subgraph isomorphism problem? 
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